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Abstract 



The hyperdeteminants considered here are the simplest analogues of determi- 
nants for higher rank tensors which have been defined by Cayley, and apply 
only to tensors with an even number of indices. We have shown in a previous 
article that the calculation of certain multidimensional integrals could be re- 
duced to the evaluation of hyperdeterminants of Hankel type. Here, we carry 
out this computation by purely algebraic means in the cases of Selberg's and 
Aomoto's integrals. 
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1 Introduction 

Selberg's integral and its generalizations are multiple integrals encountered 
in many fields such as mathematical statistics, random matrices, statistical 
mechanics, special functions, integrable systems, and can be interpreted as 
generalized moments, correlation functions or partition functions. 

Such integrals appear naturally in the computation of the normalization 
constants for various random matrix ensembles. For example the normal- 
ization of the eigenvalue probability density function for Gaussian, circular, 
Laguerre, Jacobi and Cauchy ensembles are limit cases of Selberg integrals 
(see jSH EH US] ) • Generalizations appear in the calculation of the partition 
functions of log-potential Coulomb systems [TE] . Such multiple moment in- 
tegrals are also connected with the problem of expanding even powers of the 
Vandermonde determinant 

A(ar) = det^r 1 ) = JJ (%j ~ a*) (1) 

l<i<ji<n 

in various bases of symmetric functions. Expansion of these powers in terms 
of Schur functions is still an open problem (see [H21 EH1 EI]), whose solution 
would find important applications to calculations related to some aspects of 
the fractional quantum Hall effect (see [3211101) described by Laughlin's wave 
function [2H] 




The same expansion would allow one to compute many more interesting 
integrals (see [TJ CHI 123 I2E] ) ■ 

Selberg's integral is a generalization of Euler's Beta integral, introduced 
by himself in 1944 (see [33J), in view of proving a generalization of a theorem 
of Gelfond. Selberg found that his integral could be written as a product of 
quotients of gamma functions 

5„(a,6; 7 ) = / •••/ lA^I^TTzr 1 ^-^) 6 " 1 ^ 
Jo Jo i=1 

r(a + n )T(b + n )T((j + 1) 7 + 1) 

11 r(a + 6+(n + j-l) 7 )r(7 + l) ' 1 ] 

This identity was subsequently proved in many different ways by other au- 
thors, such as Anderson PJ, Aomoto [S] and Dotsenko and Fateev ^2]- These 
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proofs are all interesting, as they involve generalizations of Selberg's integral 
or techniques which can be useful in other contexts. However, Selberg's orig- 
inal proof has the specificity that it is almost elementary, in the sense that 
he computed the integral in the case where 7 is an integer k, using only 
arguments on the degrees of the polynomials appearing after integration of 
the terms obtained by expansion of A(x) 2k , and a very simple change of 
variables. Then, he concluded by involving a classical theorem of Carlson to 
extend it to complex values of 7 (see (HI] or |HEj)- As we shall see below, the 
elementary part of this proof comes out naturally if we consider this object 
for 7 = k G N not as an integral but as a special case of a hyperdetermi- 
nant, a notion originally due to Cayley but whose only available textbook 
presentations seem to be Refs. [3IJ|35| (in Russian). 

The hyperdeterminant is a polynomial invariant of tensors of even order 
(i.e., with an even number of indices). In the case of matrices, its definition 
gives back the classical determinant. 

We have explained in [2HI that Selberg's integral could be rewritten as a 
hyperdeterminant of Hankel type (i.e., whose entries depend only on the sum 
of the indices) built from the moments of the beta distribution. In this paper, 
we give a complete translation of Selberg's proof in the hyperdeterminantal 
formalism, using only very simple tools. The hyperdeterminantal aspects of 
multiple integral computations look very promising. At present, we know 
only a few tools to handle such polynomials, but this is already sufficient 
to deal with some interesting cases. To illustrate this point, we sketch a 
hyperdeterminantal proof of Aomoto's integral 



A a n b *{y) = / •■•/ A(x) 2k T\(y-x l )xr 1 (l-x l ) b ~ 1 dx 1 ---dx n (4) 
Jo Jo i=1 

= (-2)- n S n (a,b;k)Pt l ^'\l-2y) (5) 

where P^' b are the monic Jacobi polynomials. Aomoto has proved this iden- 
tity in [3] by means of differential equations. The main interest of Aomoto's 
original proof is its relation with the Calogero model. Barsky and Carpentier 
have given in [3] a different proof, based on Anderson's method. 



2 Elementary hyperdeterminant calculus 
2.1 Definition 



We shall find it convenient to use a set of n Grassmann (that is, anticom- 
mutative, or fermionic) variables r] = {rji, . . . , rji, . . . }. A tensor of order k 
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and dimension n , that is, an element of V® k , where V is an n dimensional 
vector space, will be represented as 

M= Mi y ...i k f]i x ® • • • ® T] ik . (6) 

That is, we identify V with the vector space spanned by T]i, . . . , r\ n . The 
array (M^...^) will be called a hypermatrix. The hyperdeterminant Det(M) 
of M can be compactly defined as the coefficient of (7/1 • • -ri n )® k in 

M n = n\Det{M){r ]l ---r ]n )® k (7) 

where the power M n is evaluated in the tensor product of k copies of the 
Grassmann algebra over the n distinct ?7i's, so that the only surviving terms 
are those involving (771 • • • r] n )® k . 

For convenience, we will use the notation 

Det fe (M il ... ifc )5 l = Det(M). (8) 

Note first that if M is a tensor of odd order k, one has Det(M) = 0. If M is 
a tensor of even order 2k and dimension n, one has 

1 n 
Det(M) = - <^)---<^k)l[M (Tl(i) ... a2k{i) . (9) 

"i,->»2te6„ i=l 

Setting k — 1, we recover the classical definition of the determinant for 
matrices. 

2.2 The invariance property 

Definition (JIJ) is very useful for proving properties of hyperdeterminants. Let 
us consider the natural action of the group GL(V) x2k on V® 2k 

(g {1) ® • • • ® g m ) ■ vi ® • • • ® v 2fe = (<7« Vl ) ® • • • ® (<? (2fc) v 2fc ) (10) 

Then, for a general M e \/® 2fc , 

(^)®...®^)). M = £ .\/,,.. ; , x 

l<il,- J2fc<™ 
n n 

en) 

il=l «2fc=l 

Hence, 

Det ((gW ® • • • ® ^ 2fc >) • M) = det • • • det (g (2k ^ Det(M). (12) 



Hyperdeterminantal calculations of Selberg and Aomotos's integrals 5 



2.3 Minor summation formula 

Consider two tensors 

M= M n...i. k Vh®---®mu k (13) 

l,<i 1 ...,i 2 k<n 

and 

l<h,.. .,«2fc<™ 

Remarking that MN — NM = (in the 2/cth tensor power of the Grassmann 
algebra), one has 

n 

(M + N) n = ^2 M*N n ~\ (15) 

If I = (Ii, ■ ■ ■ , 12k) is a 2fc-tuple of subsets of {1, . . . , n}, we will denote by 
M[J] the tensor 

M[/]= M n-, 2fc %®---®^ fc (16) 

and by rji the product 

where the symbol Yli denotes the product taken in increasing order of the 
subscripts. 

Let C r nk be the set of pairs (I, J) of 2/c-tuples / = (ii, • • • ,hk) and 
J = ( Ji, ■ ■ ■ , J2k) such that for each s G {1, . . . , 2/c}, (J s , J s ) is a partition 
of {1, . . . , n} into two blocks of sizes r and n — r. For each (/, J) G fc 
we will denote by e(I, J) the product of the signs of the permutations cr s = 
(ii, . . . , i r ,ji, ■ ■ ■ ,j n ) (written as a word) where I s = {i\ < • ■ ■ < i r } and 

J s = {ji < ■ ■ ■ < jn-r}- 

We can now write, for each p < n, 

M p = p!^Det (M [I\)th (18) 

where the sum is over the 2fc-tuples / = (ii, • • • , 7 2 fc) of subsets of {1, ... , n} 
of cardinality p. Hence, by f|15|) . 

n 

Det(M + N) = ^ e(J, J)Det(M[/])Det(N[J]) (19) 
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Now, if we set 

M'= (20) 



Ki2,...do.k<n 



and 

n 

M " = J2 S M h,h,-,h k VH®---^Vh k , (21) 

ii=2 I<i2,...,i2fc<n 

remarking that M = M' + M" we obtain from 

Det(M) = Yl sign(/)M M2i ... i42fc Det(M[7]) (22) 

7=(ii=l<i2,i3v,«2fc<n) 

where sign(I) = (— l) l i+*2+-+ J 2fc anc [ / = ({1 ; . . . ; n}— ii, . . . , {1, . . . , n} — i 2 k)- 
Equations JTHJ and (J22J) can be found in [3 I34| 13 5 j . 



3 Hankel hypermatrices 

For a sequence / = (zi, . . . ,i2k), we will set |/| = i\ + ■ ■ • + i 2 k and r/j = 



3.1 Heine's integrals for Hankel determinants 

Let /i be a measure on the real line, and c n = J x n dfi(x) be its moments. 
Heine's integrals below are evaluated as Hankel determinants whose entries 
are the moments c n (see [Hi] ) 

J /\(x) 2 du,(xi) ■ ■ ■ du.(x n ) = n\ det(c i+i )o _1 . (23) 

It is easy to give a direct proof of this equality, but it can also be seen as a 
very special case of de Bruijn's integral 

/ det (4>i(xj)\ilJi(xj)) dx\ ■ ■ ■ dx n = Pf (Qij) 1 <j J < 2 „ 

(24) 

a<xi<--<x n <b 

where Qij = J [(j)i(x)i(jj(x) — (j)j(x)ip i (x)]dx and (cj)i(xj)\ipi(xj)) denotes the 
matrix whose ith row is ipi(xi), (/)i(x 2 ) , ipi(x 2 ) , . . . , (f>i(x n ), ipi(x n )\ and 

Pf (M) denotes the Pfaffian of the matrix M. In [27] , we have obtained a 
generalization of this identity, involving 2kn functions on the left hand side 
and a hyperpfaffian (see on the right hand side. In the case where the 
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determinant is an even power of the Vandermonde determinant, we obtain a 
generalization of Heine's theorem turning multiple integrals into Hankel hy- 
per determinants. This generalization is discussed in [2Hj and will be recalled 
in the following subsection. 

3.2 Generalization: Hankel hyperdeterminants 

We consider here the integral 

K,k = J ■■■ J A (^i ? • • • > x n ) 2k d^{ Xl ) ■ ■ ■ dy.{x n ) (25) 

Expanding the even power of the Vandermonde determinant, this integral can 
be expressed as a Hankel hyperdeterminant whose entries are the moments 
c n = J x n d/u,(x): 

2^ = n!Det(q fl )^ 1 . (26) 
Consider now the following two tensors 

S n (a,M)= B(a + |/|,6)77 7 (27) 

0<ii,t2,...,i2fc<n— 1 



and 



m, (28) 



0<!i,!2,...,i2fc<H-l 11/ 



where B(a, b) denotes the Beta function B(a, b) = f x a ~ x (l — x) b ~ l dx. Re- 
marking that the B(a + n, b) are the moments of the Beta distribution, one 
obtains after applying (j2T?j) 

S n (a,b;k) = n!Det (S n (a, 6; £;)) (29) 

and in the same way 

A a ^ k {y) = n!Det (A^ fc (y)) . (30) 

3.3 Expansion of Hankel hyperdeterminants 

In this section we consider the hyperdeterminant of a general Hankel tensor 

7~Ln,k = Det 2 A; (Xj 1+ ... +i2fc )p 1 . (31) 
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Expanding this as a polynomial in the X i} one finds 

«M = ( 32 ) 

A 

where the sum is over the n-tuples A = (Ai, • ■ ■ , A n ) such that |A| = kn(n— 1) 
and such that for each % e {1, • • • , n}, (i — l)k < Aj < kin + % — 2). 

This expansion is derived in [2H]- A special case appears in Selberg's 
proof (see jSUESj)- 

Let V n> k be the coefficient of X^,^ in (|3~2j): 

= c (fe(n-l))"- (33) 

If we set Xi = for i > — 1), we obtain from ()32j) that 

Det 2 fe (Xj 1+ ... + j 2 Jg 1 = t>n,fc-^fc(n-l) ■ (34) 

A closed form for D n> fc follows from a well-known constant term identity (the 
"Dyson conjecture") 1 

where C.T. means "constant term". On the other hand, 

n s \ k n i 

n C 1 -^ 1 ) "(-^rilTpny < 36 > 

i,i=i ^ 3 ' i=l X i 

Expanding the power of the Vandermonde determinant, one finds that for 



x <T 1 (i)+-+a 2k (i)-2k-k(n-l) 



C n , k = (-lfC.T.j Yl E W" ,e WlI 

= (-l) fe n!Det 2fc (5 n+ ... +l2fc , n )^ 1 (37) 

where 5ij = 1 if % = j and otherwise (Kronecker symbol). Using (j34j) . we 
obtain 

C n , fc = (-l) fc n!D n , fc . (38) 

Hence, 

fM = (-l)^f, ^ (39) 



n! V k, ■ ■ ■ , k 



1 This was conjectured by Dyson ^3 i n 1962. I. J. Good gave, in 1970, an elegant 
elementary proof involving only Lagrange interpolation |18| . 
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3.4 Minors of a Hankel hyperdeterminant 

We consider here a Hankel hyperdeterminant 

M = J2 M \i\Vh ® • • • ® ?7i 2fc (40) 

and a family of 2/c subsets J = (J 1 , . . . , J 2fc ) of a fixed size m < n with 
^ = {io ^ A ^ "' ^ im) c {0, — 1}. We want to expand the 

polynomial M[J] m , where 

m[j] = M fi % ® • • • ® ( 41 ) 

I 

One has 

Det(M[J]) = - £ ^i)- C (^)n^ l(0+ -^co- (42) 

f li-ftESm i=l 

But for each 2/c-tuple (a"i, . . . , a 2 k), one has 

II M 4(0 + -+^2 fc (0 = II M «(0+-+^ fc (0-2fc+ri l(0 +-+r» fc(0 ( 43 ) 



Z=l 



with < rj < n — m. Hence, by a computation similar to that of a complete 
Hankel hyperdeterminant, 

m m 

i=i i=i 

with — 1) < Aj < k(m + i — 2) and < Sj < 2k(n — m). And 

m 

Det(M[J])=^c^ ;J n% (45) 

A i=l 

where k(i — 1) < Aj < A;(2n — m + i — 2). 



4 Selberg's integral 

4.1 From the Hankel form to the symmetric form 

Selberg's integral admits another hyperdeterminantal representation which 
can be obtained directly from the previous one (without manipulating the 
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9 = 

l<i,j<n 



integral). It suffices to remark the following property 

n 

5^(-l) 1 ( n )B(a + z,6) =B(a,6 + ?i). (46) 
Hence, we have from (refeqlO) 

k 

S n (a, b- k) = (I, • ■ ■ , I, g~^~g) ■ S^ m (a, b- k) (47) 
where I is the identity matrix, g is the n x n matrix 

J2 (-lyf')**®^ (48) 

and 

S^ m (a,b;k) = B(o+|/|,6+|J|)^®^ (49) 

/,Je{0,...,n-l} fc 

It follows from (JT21) that 

det(#) fc Det (S^ m (a, 6; k)) = Det (S n (a, 6; k)) (50) 

and, since 

det ( 5 ) = (-I)t, (51) 

we find 

Det (S^ m (a, 6; fc)) = (-lj^^Det (S n (a, 6; fc)) (52) 
without the help of the integral representation. 

4.2 Selberg's original proof 

Each step of Selberg's proof can be viewed as a manipulation on a hyperde- 
terminant (see jHI] for example). It can be divided into two parts. In the 
first part, Selberg proved that 

c r u u\ TT T(q + jfc)r(& + jfc) 

S n (a,b k) = c nik I I — — — — - — — — . 53 
f = o r ( a + b+ + j - l)k) 

where c n ^ does not depend on a and b. In this proof, Selberg started from the 
expansion of the Vandermonde determinant. In term of hyper determinants, 
it is equivalent to the expansion (j2HJ). For our purpose, equation (jSHJ) reads 

„ , , ttt r(a + jk)T(b + jk) 

Dt(s.(».*;*)) = ^n r(a + t + ( n + j-i)t) (54) 
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where a n ^ = n\c n ^. To obtain it, we start from formula 

Det(S n (a, b; k)) = ^ c"' fc B(a + A 1; b) ■ ■ ■ B(a + A n , b) (55) 

A 

where the sum is over all n-tuples A = (Ai, • • • , A n ) such that |A| = kn(n — 1) 
and for each % G {1, • • • , n}, (i — l)k < Aj < kin + % — 2). Then, extracting 
common factors, we arrive at an expression of the form 

DetfS (a b k)) - Q{a ' b) TT T ( a + J k ) T ( b + J k ) (56) 
Det{b n {a,b,k))- R{b) Ur(a + 6 + (n + j-l)Jfe) ( } 

where Q(a,b) is a polynomial of degree at most kn(n — l)/2 and R(b) is a 
polynomial of degree kn(n — l)/2. But (JH2"j) implies that 

Det(S„(a, b; k)) = Det(S„(6, a; k)) (57) 

and then 

Q(a,b) = Q(b,a) 
R{b) R{a) 1 ' 

which implies that a n ^ = is independent of a and b. 

The second part of Selberg's proof consists in establishing a recurrence 
relation for a n 

a n ,k = r-oc n -i,k- (59) 

nl 

This can be found from the limit 



Hence, one has 



Lb-,k,n= Ihn aBet(S n (a,b;k)) (60) 



L 6;fc , n = lim Det(T a ) (61) 



where T a is the tensor defined by 



T a = aB(a+\I\,b)r] I + 

le{0}x{0,...,n-l} 2k - 1 

+ Yl B(a + I 7 b) m (62) 

,n-l}x{0,...,n-l} 2k - 1 



Remarking that 



lim a B(a + n, b) — <5 nj0 (63) 
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where is the Kronecker symbol, one obtains 

L b . k>n = Det(T) (64) 

where T is the tensor 

T = Vo ®...® m + B(|J|,6)77j (65) 

Ze{l,- ,n-l}x{0,...,n-l} 2fc - 1 



Hence, expanding by (|22|) with respect to the first subscript, one has 

L b -k,n = Det(S„_i(l, b; k)). (66) 

Together with (|54jl this equality gives the recurrence relation (|59jl and proves 
Selberg's identity. 

4.3 An alternative ending for Selberg's proof 

We shall now give a simpler proof of Selberg's identity by reducing it to the 
Dyson conjecture, which, as we have already seen, is by now a familiar and 
elementary statement involving no more than Lagrange interpolation. Using 
the functional equation of the V function and the linearity properties of the 
hyperdeterminant, we can recast the Hankel hyperdeterminantal expression 
of Selberg's integral in the form 

/ \ \ n— 1 

0)]/] 



Det (S n (a, b; k)) = B(a, b) n Det 2k v /'' . (67) 

where (a)„ = a(a + 1) • • • (a + n — 1) denotes the Pochhammer symbol. The 
result will follow if we can obtain a closed form for the hyperdeterminant 

V n>k (a,b) = Bet 2k [^fj ■ (68) 
We start as in the previous proof, obtaining an analog of ()54jl 

T>„Ja, 6) = <, f[ '")'<■-;»(*-»'<--» (69) 

„Ji {0)k{n+m-l) 

where the constant a' n k is independent of a and b. Now, for every a' 



r> I V\ TT i. a )k{m-l)ib — a) fc ( m -l) /, s.n-1 ,_„s 

V k}n {a,b) =11 — ^ i — — ^Det 2fe {(a)\i\) . 70 

(&)fe(n+m-l)(a')fc(m-l) 1/0 
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This means that Selberg's integral can be deduced from the knowledge of 
Det 2 fc((a / )|/|)Q~ 1 for any particular choice of a'. If we set a' = —k(n — 1), one 
has by flMJ) 

Det 2fe {{-k{n - l)) m )o -1 = (-l) kn ^-% n , k (k(n - 1))!™. (71) 

But d Ut k has already been evaluated by means of the Dyson conjecture (|3U|) 
and the closed form of Selberg's integral follows. 

5 Aomoto's integral 

5.1 Another hyperdeterminantal representation 

On the integral representation (@J), we obtain easily the following identity 
through the substitution 

A a f> k (y) = (-l) n A b n a -> k (l-y). (72) 

There is also a simple hyperdeterminantal proof of this assertion. We consider 
the tensor 

K a]k (±~y)= £ (l-y- a ^' /| |/| )B(6+|/|,6)7y J (73) 

/e{0,...,n-l} 2fc ^ ' 

and the action of 2k copies of the matrix g defined before. Using equality 
(Jinj), one has 

(g,...,g)A b ^ k (l-y) = -A^ k (y) (74) 

Hence, 

Det(A^ fc (y)) = (-l)"Det(A^ fc (l - v))- (75) 

5.2 Minors of Selberg's hyperdeterminant 

To prove Aomoto's identity, we need some preliminary results on Selberg's 
hyperdeterminant. Let us consider a family of 2k subsets I = (J 1 , . . . , I 2k ) 
with P = {i J < i{, ■ ■ ■ , < i J m -i} C {1, - • • , n} and the associated sub-tensor 
of Selberg's Hankel tensor: 

S n (a, b- k) [I] — B(a + i 1 jl +--- + if 2k , 6)?/^ ® • • • ® ^ . (76) 

o<ii,.--,i2fc<"i-i 

One can write by 

m 

Det(S„(a, b; k) [/]) = £ c n x kJ J] B(a + A,) (77) 

A i=l 
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where k(i — 1) < Aj < k(2n — m + i — 2). Hence, 

Det(S„(M; ,)[/])= fi r^^^^.l^) ^^^ 

_ ^ r(g + ib(t-i))r(6) 

where p/(a, 6) is a polynomial. In the same way, we will use the following 
identity 



n - m 

D*(S.<« + 1, 6; *)[/]) = n „ + b + k{2n ._ j _ t) + *<«-!» 



X 



n m 



x n (n + 1 + H2n n rgFrfeS)) ^ + »>• (79) 

which follows immediately from (|78j). 

5.3 A proof of Aomoto's identity 

We start from the hyperdeterminantal representation and we remark that 

A.f>\y) = S(a, b; k)y - S(a + 1, b; k). (80) 
From ()19p. we obtain 

n 

Det(A***(y)) = ^(-l)"" r ^ x 

r=0 

x ^ e(J, J)Det(S(o,6; fc)[/])Det(S(o + 1,6; fc)[J]). (81) 
From (jZSj) and (JZHJ), one has 

Det(A^)) = n a + 6 + fc( L-^-l) X 

v „ r n r(o + fc(t-i))r(6) 

^ llr(a + 6 + ife-2n-i-l) X (§2) 



n n 

X ' 

r=0 i 

n—r 

X 

j=l i=n— r+1 



J[(a + fc(i-l)) ] [ (a + 6 + &(2n - i - l))<p r (a, 6) 
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where ^3 r (a, b) is a polynomial in a and b. 

On the other hand, we expand the hyperdeterminantal representation of 
Aomoto's integral by (J32|) and obtain, after extracting common factors as 
above, 



Det(A^ fc ) = Bet(S n (a,b ] k))J2^U{y^-i ) ( a + b ^ 



n,k n 

k(n+i-2)-l 



+ 6)2S*" 30 " 1 } (83) 



where Q(b) is a polynomial of degree kn(n — 1) in b, (a)™ denotes (a + n)(a + 
n + 1) • - • (a + m) if n < m and (a)™ = 1 otherwise and the sum is over the 
partitions A = (Ai < • • • < A n ) verifying kii — 1) < A» < k(n + i — 1) for each 
j G {1, • • ■ , n}. After combining the coefficients of the y l7 s, one obtains 

Det(A^(y)) =Det(S n (a,6;fc))^2/ r ^^ (84) 



r=0 



where r (a, 6) is a polynomial. Hence, comparing the coefficients of y l in the 
expressions (|H2J) and (|H1|) . we find 

n— 1 1 

Det(A^*(y)) = Det(S n (a, 6; fc)) JJ - - - - ^ - . _ ^ x (85) 



i=0 



x E^-^) 1 n( a + ^-^) II (a + b + k(2n-j-l)) 



i=0 tv,/ j=i j =n _j_|_i 



where Pj(a, 6) is a polynomial. Remark that each Pi(a, b) is of degree in b at 
most kn{n — 1). Now, we apply (|75|) and we equate the coefficients of in 
the left and right hand sides. After simplification, one obtains for each i 

itfy PAb,a)F 3 (b,a)Q(b) = (-lf^a, b)ft{a, b)Q(a)) (86) 



where 



Fi{a,b) = Y[{a + k{n - m)) J| (a + b + k(2n - m - 1)). (87) 

m=l m=n— i+1 

Remarking that does not divide Fi(a,b), equality (jHB|) implies 

Pi{a,b) = Ri{a)Q{b). (88) 
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Hence, setting as in the previous section, a = — k{n — 1) in 

(j) #i(W(M = (-l) n ' l R t (a)F t (a,b) (89) 
and remarking that 

Fi(-k(n- l),b) = (90) 

and 

n 

Fi{b,-k{n-l)) = Y[(b + k{n-m)), (91) 



m=l 

one obtains the recurrence relation 

n 



R l (b) = - J2 (-)^( & )- 



(92) 

3=1+1 v ' 

The starting point of the induction is R n {b) = -j, which follows from the 
fact that the coefficient of y n in the expansion of the hyperdeterminant is the 
value of Selberg's integral. Solving the recurrence relation, one finds that 
each Riib) is in fact independent of b, and that 

m = tgi r\ (93) 



1 

x 



Hence, 

n— 1 

Det(A^ fe (y)) = Det(S n (a, 6; fc))— 



n! -!- a + 6 + kin + i — 1) 

1=0 

Y[(a + k(n - m)) JJ (a + 6 + k{2n - m - l))y\M) 



i=0 

n n—i n 

X 

i=0 m=l m=n— i+1 



. 1 )»Det(S„ (a , 6; t)) ? ^ , ± + i + » " ^ 

V fc k In i=0 



We recognize a hypergeometric function of type iF\ which can be evaluated 
as a monic Jacobi polynomial (see [2] for details) 



Det(A^*(!/)) = (-irDet(S B (a,6;fc)) ^ ^ x 



x 2 ^i 



-n,n + f + |-1 



(95) 



n! 

Aomoto's equality is therefore proved. 



-2")"™ a_i A_i 

-Det(S„(a, (1 - 2j/) 
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6 Conclusion 

The examples discussed in this article show that the hyperdeterminantal 
calculus is a pertinent tool to handle an interesting class of multiple integrals. 
The hyperdeterminant is a particular invariant of hypermatrices. In the case 
of antisymmetric hypermatrices, another invariant has similar properties: 
the hyperpfaffian &T\. In the classical case of matrices, de Bruijn has 
shown in [TTj that multiple integrals of some determinants could be evaluated 
as Pfaffians. In [2Zj, we have found some generalizations to hyperpfaffians 
of these identities. In fact, the generalized Heine theorem is a particular 
case of the generalized de Bruijn integral, where the determinant can be 
factored into a product of determinants of matrices whose dimension is the 
number of integration variables. There exist other invariant polynomials 
of hypermatrices |HJ El HIH Ej which can be, in principle, computed using 
Cayley's Omega process or other methods of invariant theory (see U\ 
for examples involving 2x2x2x2 and 3x3x3 hypermatrices) . A natural 
question is whether there exist other integral identities involving some of 
them. 

Other generalizations of Selberg integral are encountered in the physical 
litterature. For example the Dotsenko-Fateev (see ^2]) and Kaneko integrals 
(see [2PJ give the partition functions of various systems. 

Kaneko's integral reads 

/in/ n \ 

|A(x)| 2 ^n Xi^Jlfa - Vj )dxi J (96) 

and is evaluated in terms of generalized orthogonal polynomials associated 
with Jack polynomials (see [2H1 EU EH] ) ■ By the generalized Heine theorem, 
when 7 is an integer, Kaneko's integral can be evaluated as a hyperdetermi- 
nant of moments (see [28J). Kaneko's proof of his identity is related to the 
Calogero- Sutherland model. Up to now, we have been unable unable to find 
a hyperdeterminantal interpretation of this. In the same way, Anderson's 
(see [l]) and Aomoto's (see jH|) proofs of Selberg's integral seem to contain 
information unrelated to the hyperdeterminantal representations. 
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